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Abstract. - The dewetting dynamics of ultrathin polymer films, e.g. in the model system of 
polystyrene on a polydimethylsiloxane-covered substrate, exhibits interesting behavior like a fast 
decay of the dewetting velocity and a maximum in the width of the built-up rim in the course 
of time. These features have been recently ascribed to the relaxation of residual stresses in 
the film that stem from the nonequilibrium preparation of the samples. Recent experiments by 
Coppee et al. on PS with low molecular weight, where such stresses could not be evidenced, 
showed however similar behavior. By scaling arguments and numerical solution of a thin film 
viscoelastic model we show that the maximum in the width of the rim can be caused by a temporal 
evolution of the friction coefficient (or equivalently of the slip length), for which we discuss two 
possible mechanisms. In addition, the maximum in the width is affected by the sample age. As 
a consequence, knowing the temporal behavior of friction (or slip length) in principle allows to 
measure the aging dynamics of a polymer-polymer interface by simple dewetting experiments. 
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Introduction. — Dewetting experiments [1-3] have 
been proven to be a simple tool to investigate thin poly- 
mer films on various substrates [4,5]. Especially the model 
system of the incompatible polymers polystyrene (PS) 
on polydimethylsiloxane (PDMS) grafted on a silicium 
substrate has attracted numerous experimental investiga- 
tions. There is consensus that friction and slippage at the 
interface, as well as the viscoelasticity of the dewetting 
film are crucial to understand the experiments [6,7]. To 
account for the occurrence of a maximum in the width 
of the rim in the course of time, and to interpret the ef- 
fects of sample aging on this peculiar feature of the thin 
films' dewetting dynamics, residual stresses in the PS film 
have been shown to be important [8-10]. These stresses 
are supposed to stem from the sample preparation process 
(usually spin-coating) and represent an additional degree 
of freedom that can (partially) relax during dewetting. For 
high enough molecular weights, residual stresses have been 
clearly evidenced in various systems [8,11,12]. However, 
their relaxation well below the glass temperature, needed 
to interpret the effects of aging [8], is currently debated. 
An additional puzzle is the fact that for films of low molec- 
ular weight no such relaxation of stress could be evidenced, 
while a maximum in the rim width and its dependence 
on aging have still been found. Very recently, in parallel 



to dewetting experiments the interface between the two 
polymer species has been studied by neutron reflectome- 
try [13] and interdiffusion of the PS-PDMS interface has 
been found, for samples aged well below the PS bulk glass 
temperature. 

In this letter we show by scaling arguments and demon- 
strate by numerical solution of a thin film model that a 
slow increase in the friction coefficient or, equivalently, a 
decrease of the slip length can give rise to a maximum in 
the rim width. This evolution of friction/slip might be due 
to several mechanisms: first, as suggested by the neutron 
reflectometry measurements [13], a slow roughening of the 
film-substrate (PS-PDMS) interface with concomitant in- 
crease in friction may be responsible. A second possibility 
would be that (for the lowest PDMS grafting densities) a 
small number of PS chains may enter the PDMS brush and 
become attached to the substrate. It is known that minute 
amounts of such 'connectors' [14-16] already decrease the 
slip length. 

A second finding is that the value of the maximum 
rim width decreases monotonously with the friction co- 
efficient at the beginning of the dewetting process. This 
could explain the experimentally observed decrease of the 
maximum with the sample age. Thus knowing the evo- 
lution of the friction coefficient allows in principle to re- 
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late the maximum width to the aging history. As stud- 
ied in Refs. [8,9,17,18], a similar connection between rim 
width and aging time can also be caused by relaxing in- 
ternal residual stresses. One would expect that stress re- 
laxation is the dominating mechanism for high molecular 
weight films, while for low molecular weight films interface 
changes are responsible. In general, the two nontrivial pro- 
cesses, stress relaxation in bulk and evolution of the inter- 
face, might be present at a time, which leaves interpreta- 
tion of experimental results an intricate task. The specific 
dynamics of the PDMS brush-PS melt interface being un- 
known, drawing conclusions from the dependence of the 
rim width on the sample aging, as suggested above, is 
only a rough estimate at present. Dynamic self-consistent 
field techniques or molecular simulations would be needed 
to establish this dynamics. 

Which processes can lead to a maximum in the 
rim width ? A scaling argument. First we want 
to give a simple argument, why a slow evolution of the 
interface makes the occurrence of a maximum in the rim 
width possible. For dewetting of a viscoelastic film, scal- 
ing arguments have been used to qualitatively understand 
the dynamics [17,19]. The balance of the work done by 
the driving force and the dissipation by friction reads ap- 
proximately \S\V(t) ~ CW{t)V(tf or 

l|l ~W(t)V(f). (1) 

Here V(t) is the velocity at the dewetting edge and W{t) is 
the width of the rim. The driving force acting on the rim 
is the negative of the spreading parameter S = 7™ — 7^ — 7 
(where the 7's are the interface energies for the substrate- 
vapor, substrate-liquid and liquid-vapor interfaces, respec- 
tively), and C is the friction coefficient associated to the 
polymer-polymer interface. In general, the friction will be 
nonlinear (see below), but for the simple argument given 
here this is not of importance. In addition to this bal- 
ance, one needs the fact that the dewetting velocity is 
monotonously decreasing - there is no mechanism, and also 
no experimental evidence, for a speed-up of the dewetting 
process. 

What possibilities does Eq. (1) allow for W(t) to have a 
maximum? Clearly, if W has passed through a maximum, 
it has to decrease for some time simultaneously with the 
monotonously decreasing velocity. This is not possible for 
Eq. (1) if the left hand side is a constant. A necessary 
condition for a maximum in the width is thus that the 
effective driving force \S\/( decreases in the course of time. 
There are several possibilities to achieve such a decrease: 
The first one is motivated by residual stresses stemming 
from spin-coating. Introduction of a residual stress a(t), 
that is allowed to relax in the course of time, leads to 
a renormalization of the driving force IS"! — > \S\ + hoo~(t) 
[17,18] (where h is the film thickness). Accordingly upon 
stress relaxation the driving force decreases. 

The second possibility, suggested by the interface evolu- 



tion evidenced in Ref. [13], is an increase in the friction co- 
efficient C — ► £(i) due to roughening of the liquid-substrate 
interface. Additional to the increase in friction, interface 
roughening also leads to a decrease of the driving force 151, 
as discussed for autophobic dewetting in [20] , and can even 
lead to a cross-over to wetting in that case. For the incom- 
patible system PS-PDMS, where the PS-air surface energy 
is dominating, this latter effect is however not probable for 
the expected microscopically small roughnesses. 

A third possibility, also associated to tiny changes in 
the interface, is the successive attachment of a few PS 
chains to the substrate where the PDMS-coverage is low 
or defective. This process leads to a decrease in the driving 
force IS"! — — i/(t)lf, where v is the areal density of such 
'connector' chains, I the length they are stretched upon 
pull-out and / the pull-out force [14, 15]. Additionally 
the friction increases due to the connectors (in simplest 
approximation £ — > ^ + v(t)K with k the friction coefficient 
of connectors; for a more accurate treatment, sec [21]). 

In the following we would like to exemplify the second 
case further. We first give an estimate of how the interface 
evolution gives rise to an increasing friction coefficient. 
Later on this is implemented in a thin film model and its 
influence on the dewetting is studied. 

Interface evolution. — Polymer-polymer interfaces 
are well studied [22,23], although mostly at equilibrium. 
The equilibrium width of an immiscible interface is given 
by S oc a/^/x [24], where a is the Kuhn length and \ 
is the Flory interaction parameter. A simple scaling ar- 
gument [25], consistent also with a full statistical treat- 
ment [26], states that the slip length b at the interface 
is given by the ratio of bulk and interfacial viscosity and 
scales as b oc N/5 2 . The factor N is the degree of poly- 
merization and leads to the well known large slip lengths 
for polymers [27]. Alternatively, the slippage length can 
be defined as b = implying that the friction coeffi- 
cient scales like ( oc 5 2 . This scaling should hold for the 
interface between two immiscible melts, for a brush-melt 
interface the situation is slightly different. The structure 
of brushes in various matrix polymers has been studied in 
Ref. [28], and the width for the immiscible case has been 
found to be narrow, but larger than the width between 
melts. The effect of the grafting density on the interface 
width is unclear, but one would expect that the grafting 
density of the PDMS-brush should not be too high to pre- 
vent total stretching and allow for mobility of the brush 
chains. Indeed, as reported in Ref. [13], roughening was 
pronounced only for the lowest grafting densities where 
the interface width increased from 1 to 4 nm in the course 
of time. It should also be mentioned that polymer in- 
terfaces display thermally activated capillary waves [29], 
which lead to a higher apparent interface width. 

As far as the dynamics is concerned, the evolution of 
the interface towards its equilibrium value due to chain 
interdiffusion and its influence on the dewetting via the 
slip length or friction is difficult to determine. Strictly 
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speaking it requires to solve a hydrodynamic problem for 
a system with an evolving boundary, to be determined by 
a model for an immiscible polymer brush-polymer melt 
interface and which in turn will be influenced by the flow. 
This scope clearly demands extensive microscopic mod- 
cling and we rather use a simple estimate, motivated by 
experiments. Intcrdiffusion dynamics for the symmetric 
case has been studied for both a melt and a network of 
PS on a PS-brush. The interface width has been found 
to change by a factor of 3-4 on the time scale of a repta- 
tion time [30,31] and displayed a logarithmic behavior in 
time. For the immiscible system PS-PMMA, also a loga- 
rithmic growth in time has been found [32], and has been 
explained by capillary wave theory. For the purpose of 
demonstrating the influence on the dewetting, we will as- 
sume a logarithmic time dependence of the roughness. In 
combination with C oc S 2 , as a rough estimate for the time 
dependence of the friction coefficient we thus suggest the 
following expression: 



C(*) = Co + Ci iog 2 (i + t/n 



(2) 



Here Co is the friction coefficient associated to the status 
of the interface upon start of the dewetting process (by 
heating above the PS-glass temperature), £i governs the 
importance of the temporal evolution and t* is a short 
time scale (the time t* has to be introduced on dimen- 
sional grounds; we assume t* <C To, t± as defined in Eq. (4) 
below). We admit that the dependence given by Eq. (2) 
is empirical and a microscopic foundation would be desir- 
able. However, as shown below, the effects of an increasing 
friction on dewetting turn out to be very robust with re- 
spect to variations of the detailed form of 

Modeling. — In order to demonstrate that a time 
evolution of interfacial properties indeed can give rise to 
a maximum in the rim width and a concomitant depen- 
dence on sample age, we performed a numerical study for 
the model established recently in Rcfs. [9,17,18]. This 
model is a lubrication approximation for a Jeffrey-type 
viscoclastic fluid film on slippery substrate. As motivated 
by recent experiments on polymer-polymer friction [33], 
nonlinear friction with the substrate is considered. In the 
model, one accounts for three fields: the velocity in the 
film v(x,t), the stress field a(x,i), and the height profile 
h(x,t) [assumed to be small compared to the slip length]. 
A simple edge geometry is considered, where the edge of 
the film is initially at x = and the dewetting occurs in 
+x-direction. As governing equations one has the momen- 
tum equation (3), a Jeffrey- type constitutive law (4), and 
volume conservation (5) 



7 + T\dt<J 
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(3) 
(4) 
(5) 



Eq. (3) is the balance of the frictional force at the film- 
substrate interface and the divergence of the total stress 



inside the film. The exponent a € [0, 1) characterizes the 
nonlinear behavior of friction and is chosen to be a = 0.8, 
a value obtained from recent experiments on high molec- 
ular weight PS dewetting on PDMS-covered substrates [9] 
and consistent with measurements of rubber-brush fric- 
tion [33]. v is a characteristic velocity that enters the 
characteristic scales, see Eq. (6) below. In Eq. (4), G is 
the elastic modulus and To and n, with to *C ti, are two 
characteristic time scales (with t\ of the order of the rep- 
tation time). These define a short-time viscosity 770 = Gto 
and a long-time viscosity 771 = Gt\ 3> As boundary 
conditions, at the edge of the film the height-integrated 
stress has to equal the driving force, h(L)a(L) = —\S\. 
Additionally we impose v(x=oo)=0, i.e. the film is unper- 
turbed far away from the edge, and as initial conditions a 
quiescent film of thickness hg. 

Solving Eqs. (3-5) analytically for short times, one can 
establish the characteristic velocity scale Vb, Q and the 
characteristic width of the rim Wq a 



V 2 



(6) 



Here c = {^-^) 2 " is a numerical prefactor and 



Vn = 



I s\ 



Wo = 




(7) 



are the scales for a linear friction law; in that simple 
case, the characteristic width can be rewritten using the 
slippage length, Wo = Vhob, an expression already pro- 
posed some time ago [34] . A rescaled version of the model 
(x 1 = x/Wo, a , v' = v/Vo, a , f = t/ r o) can then be solved 
by a shooting method on a moving grid in space and Eulcr 
iteration in time. For more details on the model and its so- 
lution we refer to [17,18]. There, the model has been stud- 
ied in detail in the presence of a homogeneous initial resid- 
ual stress (i.e. with initial condition a(x, t=0) = do ^ 0), 
that relaxes upon dewetting and thus giving rise to a maxi- 
mum in the width. Assuming low molecular weight for the 
PS film, here we neglect the possibility of residual stress 
(a(x,t=0) = 0). Instead, motivated by the roughening of 
the interface found in [13] that should lead to increasing 
friction with the substrate, in Eq. (3) wc implemented the 
time-dependent friction coefficient ((t) as given by Eq. (2). 

To test the robustness of the dewetting behavior upon 
variations in the interface dynamics, in addition we inves- 
tigated a friction evolution like 



C(*) = Co + (Coo - Co) [1 - exp(-i/T*)] 



(8) 



This dependence is motivated by the fact that the rough- 
ness, and thus the friction coefficient, should saturate for 
long times, both because of the enthalpy cost given by \ 
and due to the grafting of the PDMS chains. 

Numerical results for the temporal evolution of the 
width of the rim are shown in Fig. 1. One can see that 
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Fig. 1: The width W of the dewetting rim as a function of 
time in semi- logarithmic scale and for different initial values of 
the friction coefficient, Co- In panel (a) the friction coefficient 
evolves as given by Eq. (2) [parameters t* = 0.1, Ci = 0.03] 
and in panel (b) by Eq. (8) [parameters tq = 500, Coo = 5]. 
Viscoelastic time scales To = 1, Ti = 100 and friction exponent 
a = 0.8. Co = 1 (black), 1.5 (red), 2 (green), 2.5 (blue). The 
dashed curves in panel (b) have been obtained with constant 
friction coefficients C — Co = 1 (upper curve) and C = Coo = 5 
(lower curve). 



Fig. 2: (a) The velocity V at the edge in double logarithmic 
scale for the friction evolution given by Eq. (8) [colors cor- 
respond to the initial values of Co as used in Fig. 1(b); same 
parameters]. A rapid initial decay of the velocity is found, only 
weakly depending on Co- (b) The maximum of W , W* , as a 
function of Co for the friction laws given by Eq. (2) [Ci = 0.03 
(□) , Ci = 0.06 (0); t* = 0.1] and by Eq. (8) [r* = 500 (o), 
r* = 200 (a); Coo = 5]. Solid lines are fits to the estimate given 
by Eq. (9). 



the two different evolution laws for the friction coefficient, 
Eqs. (2) and (8), amount only to quantitative differences. 
The mechanism is robust and the qualitative behavior, (i) 
the occurrence of a maximum in the width and (ii) its 
monotonous decrease as a function of increasing Co, is un- 
changed. For the friction law with saturation, Eq. (8), it 
is instructive to compare with the cases of non-evolving 
friction coefficients fixed at Co an d Coo, see the upper 
and lower dashed curves, respectively, in Fig. 1(b). In 
both cases, there is no maximum in the width. The two 
curves enframe the ones with evolving friction. The curves 
with evolving friction will finally join for long times with 
the curve obtained with Coo (however, the present model 
should not be applied for such long times). Fig. 2(a) shows 
the dewetting velocity as a function of time. A fast decay 
of the velocity in the elastic regime, for tq < t < ti, is 
obtained as also observed experimentally. This signature 
is due to the decrease of the 'effective' driving force \S\/(, 
and has also been observed in the model when including 
a relaxing residual stress [17, 18]. 

Let us now discuss the effect of aging. Even if kept well 
below the bulk-PS glass transition, the PS-PDMS inter- 
face evolves towards its equilibrium width, as evidenced 
in [13]. This is most probably due to the fact that motion 
at the interface is easier than in bulk and that, for the 
aging temperatures used, the PDMS is not glassy. Con- 
sequently this means that aging time translates into in- 
creasing values of Co, i- e - of the value of the friction coeffi- 
cient at the moment where the dewetting process is initi- 
ated (experimentally by increasing the temperature above 
T g ). As one can see from Fig. 1, increasing Co leads to 
a monotonous decrease in the maximum of the width, in 
accordance with the experiments [13]. Fig. 2(b) shows the 



dependence of the maximum in the width as a function 
of the initial value of the friction coefficient, Co- There is 
consistent decrease for both proposed evolutions for the 
friction coefficient and for various parameter values. 

The numerically obtained results can also be interpreted 
by a scaling estimate: the typical scale for the character- 
istic width, Eqs. (6,7), implies a dependence on viscosity 
r\ and friction C as 

W oc ry^f Q-ih; ; (9) 

where a is the nonlinear friction exponent. At short times, 
one has to use rjo and Co in this formula, while at times 
t > ti one has 771 and C(i). This provides a simple picture 
for the occurrence of the maximum: The increase of the 
viscosity leads to an increase in the width of the rim up 
to a time of the order of t\. Later on, the width will 
decrease slowly since the viscosity stays at 771, while the 
friction coefficient C still slowly increases. Assuming the 

maximum W* to be roughly at t\, ci. Fig. 1, one can 

1 

estimate W* tx [C( r i)] 2 ~° • This expression was used to 
fit the dependence of W* as a function of Co to obtain the 
solid lines in Fig. 2(b), with a fixed value a = 0.8 for the 
exponent, as used in the numerical solution. Although the 
window of Co _va lues is too small to establish a real power 
law, the agreement of this simple estimate is very good. 

The numerical study and the simple scaling thus suggest 
that - once the temporal increase in friction due to rough- 
ening is established by a microscopic theory - it would 
be possible to measure the roughness evolution below the 
glass temperature, i.e. the aging of the buried polymer- 
polymer interface, by tracing the maximum in the rim 
width in a dewetting experiment. 

We should mention that in case of the interface evo- 
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lution studied here, a maximum in the width of the rim 
can already be obtained for a linear friction law, i.e. for 
a = in the momentum equation (3). This is in contrast 
to the case where relaxing residual stresses are responsible 
for the maximum, as discussed in detail in Refs. [17,18], 
where nonlinear friction was crucial. The underlying rea- 
son for this difference is the fact that the friction coeffi- 
cient already influences W for a linear friction law, see 
Eq. (9) with a = 0, while for the case with residual 
stress the nonlincarity in the friction was needed to in- 
troduce a dependence on residual stress in the scaling for 
W . However, looking at the scaling law Eq. (9), for a = 
the dependence on the friction coefficient is much weaker, 
W oc C _1 ^ 2 i as compared to W oc (~ 5 / 6 for a = 0.8. Hence 
a much larger increase in friction is needed if the friction 
law is linear, which might be improbable for the small in- 
terface changes that are to be expected for an immiscible 
polymer-polymer interface. Thus, although nonlinear fric- 
tion is not essential it amplifies the effect of the interface 
evolution with respect to the viscosity evolution due to 
viscoelasticity, as directly shown by Eq. (9). Investigating 
the (non-) linearity of the friction experimentally, similarly 
as has been done in Ref . [9] , might thus be helpful to bet- 
ter understand quantitative features of the maximum in 
the rim width. 

Conclusions and perspective. To conclude, an 
increase in friction (or a decrease in the slip length) due 
to an evolution of the substrate-film interface can give 
rise to a maximum in the rim width upon dewetting. 
This implies that the - macroscopically observable - max- 
imum in the width may be caused by and is susceptible 
to minute changes of the interface. In addition, there is a 
monotonous relation between the maximum of the width 
and the friction coefficient at the beginning of dewetting. 
Concerning the experiments on PS-PDMS systems, this 
friction coefficient in turn is affected by the aging time of 
the polymer-polymer interface below the glass tempera- 
ture. As compared to the relaxation of residual stresses 
put forward recently [8,9,17], our work offers an alter- 
native interpretation of these aging effects, in the case of 
low-molecular weight dewetting films. However, at present 
stage it is not possible to give a direct relation between the 
maximum in the rim width and the aging time, because 
neither the exact time evolution of the interface towards its 
equilibrium value, nor the exact relation (for an asymmet- 
ric melt-brush system) between interface width/roughness 
and friction/slippage is known. We would like to encour- 
age MD simulations or dynamic self-consistent field stud- 
ies to establish this dependence. Then it would be possible 
to investigate buried polymer-polymer interfaces by such 
relatively simple and cheap experiments as dewetting. 
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